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O i We briefly outline the two popular approaches on radiative corrections to 

^ I neutrino masses and mixing angles, and then carry out a detailed numerical 

■ analysis for a consistency check between them in MSSM. We find that the 

p ■ two approaches are nearly consistent with a discrepancy of a factor of 13% 

^ . in mass eigenvalues at low energy scale but the predictions on mixing angles 

D i are almost consistent. We check the stability of the three types of neutrino 

^ ' models, i.e., hierarchical, inverted hierarchical and degenerate models, under 

radiative corrections, using both approaches, and find consistent conclusions. 
^ ■ The neutrino mass models which are found to be stable under radiative 

corrections in MSSM are the normal hierarchical model and the inverted 
hierarchical model with opposite CP parity. We also carry out numerical 
analysis on some important conjectures related to radiative corrections in 
MSSM, viz., radiative magnification of solar and atmospheric mixings in case 
of nearly degenerate model having same CP parity (MPR conjecture) and 
radiative generation of solar mass scale in exactly two-fold degenerate model 
with opposite CP parity and non-zero f/es (JM conjecture). We observe 
certain exceptions to these conjectures. Finally the effect of scale-dependent 
vacuum expectation value in neutrino mass renormalisation is discussed. 
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1 Introduction 



Recent developments in the determination of neutrino masses and mixing an- 
gles from various oscillation experiments, have strengthened our knowledge of 
neutrino physics[l]. In order to have a meaningful comparison of the theoret- 
ical predictions on neutrino masses and mixing angles within the framework 
of GUTs with or without supersymmetry, with the datas from various neu- 
trino oscillation experiments [2], the effects of radiative corrections are very 
essential [3]. Considerable progress has been achieved at this front and this 
can be mainly classified into two categories - (i) the evolution of the renor- 
malisation group equations (RGEs) from high to low energy scale[4,5,6,7], 
and (ii) the low energy threshold corrections [8]. In case of running the RGEs 
in (i), the general underlying motivations are : to check the stability of the 
neutrino mass model under radiative corrections [9, 10]; to generate solar mass 
scale and also reactor angle jC/esI from radiative corrections [11]; to generate 
correct radiative magnifications of solar and atmospheric mixing angles from 
CKM-like small values at high scale[12]; and to get suitable deviations from 
the bimaximal solar and atmospheric mixings through radiative corrections 
[13,14], etc. One basic difference of the first one from the last three points 
is the very definition of the stability criteria under radiative corrections [15]. 
Radiative stability generally means that the effects of radiative corrections do 
not alter substantially the good predictions on neutrino masses and mixings 
already acquired through seesaw mechanism at high scale Mr. 

Within the framework of running the RGEs from high energy scale to 
low energy scale, we have again two different approaches so far employed in 
the literature. In the first approach (we call 'Method A' for simplicity) the 
running is carried out through the neutrino mass matrix mn as a whole, 
and at every energy scale one can extract neutrino masses and mixing angles 
through the diagonalisation of the neutrino mass matrix calculated at that 
particular energy scale[6,7,10]. In the second approach (we call 'Method B') 
the running of the RGEs is carried out in terms neutrino mass eigenvalues 
and three mixing angles [16,17,18]. We confine our analysis to CP-conserving 
case, neglecting all CP phases in the neutrino mixing matrix. 

In the present paper we carry out a detailed numerical analysis of these 
two aprroaches for a consistency check on numerical accuracy, and find out 
the stability criteria of the main three neutrino mass models[19]. We give all 
the zeroth-order as well as full textures of the left-handed neutrino mass ma- 
trices obtained from seesaw mechanism in Appendix (A,B), and use these ex- 
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pressions for checking the stabihty criteria. In addition, we further study the 
vahdity of some existing conjectures based on radiative corrections. The ef- 
fect of running vacuum expectation value to the evolution of neutrino masses, 
is further examined in both approaches. The paper is organised as follows. 
In section 2, we briefly outline the main points of the two approaches on 
renornalisation group analysis. The numerical analysis and main results are 
presented in section 3. Section 4 concludes with a summary and discussion. 

2 Renormalisation group analysis of neutrino 
masses and mixings 

We present a very brief review on the two main approaches of taking quan- 
tum radiative corrections of neutrino masses and mixings in MSSM. Our 
main motivation is to have a numerical consistency check on the results of 
these two approaches, and apply again to check the validity on some existing 
conjectures related to radiative corrections. 

2.1 Method A: Evolution of neutrino mass matrix 

In this approach the quantum radiative corrections are taken on all the el- 
ements of the neutrino mass matrix rriLL in the basis where charged lepton 
mass matrix is diagonal. The diagonalisation of the neutrino mass matrix at 
any particular energy scale leads to the physical neutrino mass eigenvalues as 
well as three mixing angles. The neutrino mass matrix mLL{t) which is gen- 
erally obtained from seesaw mechanism, is expressible in term of K{t), the 
coefficient of the dimension five neutrino mass operator in a scale-dependent 
manner, t = ln{iJ,/lGeV), 

mLL{t)^vlK{t) (1) 

where the vacuum expectation value (vev) is Vu — vq sin /3, and vq — 174GeV 
in MSSM. The evolution equation of the coefficient K{t) in the basis where 
the charged lepton mass matrix is diagonal, is given by[10] 

jK{t) = - X [^1 + Qgl - Qh] - h%, - h%,] (2) 

We replace K(t) by mLL{t) in the above equation where we assume that vev 
is scale-independent. Upon integration from high {B — L) breaking scale 
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tii{— \n{Mji/lGeV) to top-quark mass scale to{— ln{mt/lGeV)) where to < 
t < tR, we get 



rriLLito) = 



V m3i(ifl)e-^-(*o) m32(iij)e-^-(*o) m33{tR)e-^^-^'°^ J 



Ro, (3) 



Here the overall factor Rq which does not affect the mixing angles, is given 
by the expression, 

Ro = exp[{6/5)Ig,{to) + 6/^2(^0) - 6/t(to)] (4) 

which strongly depends on top-quark Yukawa coupling and gives very large 
contribution for small tan/? value (=1.60). The integrals in the above ex- 
pressions are defined as[7,10] 

1 2, 



^..(^o) = 9t{t)dt (5) 



and, 

where i = 1,2,3 and f = t,b,T respectively. The numerical values of these 
integrals at different energy scales, can be calculated from the running of 
the RGEs for three gauge couphngs (91,92,93) and the third family Yukawa 
couplings (/it, hb, hr) in MSSM. The mass eigenvalues and the MNS mixing 
matrix [20] arc estimated through diagonalisation of rriLL at every point in 
the energy scale, 

and Vmns = ^II which is identified with IJ in the neutrino oscillation 
relation 

>= U^i\vi > (8) 
where / = r, //, e and i = 1, 2, 3. The MNS mixing matrix [20] 



U 



MNS 



I f/el f/e2 f/e3 

U,, U,2 U^3 I (9) 

V Url Ur2 Ur3 
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is usually parametrised in term of three rotations (neglecting CP violating 
phases) by 



Umns — 



' C13C12 C13S12 Si3 

~C23Si2 — C12S13S23 C12C23 — S12S13S23 C13S23 I (10) 
V S12S23 — C12S13C23 — C12S23 — C23S13S12 C13C23 



where s^j = sin Qij and Qj = cos Bij respectively. The unitarity conditions are 
also satisfied 

Hence tan ^12 = |C42|/|C^ei| and tan ^23 = |C^//3|/|CA-3| and sin ^13 = |C/e3|- The 
solar LMA MSW solution favours the 'light-side' tan ^12 < 1 of the data 
[21,22,23] for the usual sign convention |m2| > It is also possible to 

express mixing angles in terms of sin^jj directly[24]. 

2.2 Method B: Evolution of neutrino mass eigenvalues 
and mixing angles 

Here we follow the main expressions from Ref.[16] for the evolution of three 
neutrino mass eigenvalues and mixing angles^. In this approach the RGEs 
for the eigenvalues of coefficient Kit) can be expressible as 

= E - + errihl) + 2hpi\K, (11) 

Neglecting /i^ and compared to h^, and taking the scale-independent vev 
as before, we have the complete RGEs for three neutrino mass eigenvalues [16], 

^ T^t^"^^' " + ^^'^ + '^^'Ul^]m^ (12) 

One can have an approximate analytical solution of the above equation by 
neglecting the small effect due to the change of C/^j in the integration range, 
as[15] 

6 

mi{tQ) = mi{tR)exp{-Igi + 6Ig2 - 6It)exp{-2U^-Ir) (13) 



^In Ref.[18] the RGEs for the neutrino masses, mixing angles and CP phases are also 

derived, whereas in Rcf.[17] a thorough discussion on fine-tuning between initial conditions 
and radiative corrections is given for quasi-degenerate neutrino masses for two generations. 
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These equations lead to the equations derived in JM conjecture [11] and have 
interesting consequences. However it should be emphasised that these equa- 
tions are valid when the mixing angles are only static. The corresponding 
evolution equations for the MNS matrix elements Ufi are given by [16] 

where f — e, ii,t and i — 1, 2, 3 respectively. Neglecting hf^ and as before, 
one can simplify the terms, 

Denoting Aui — sik+m -^^e can write the RGEs for all the elements of MNS 
matrix. For example, we give here only three of them relevant to our require- 
ment, 

^[Ur2K{A32UezUl + A^^U^lUl)] (15) 

[U^^hliA^^uMr + MzU,2Ul)] (16) 

[U.^hl^A^^U^Mr + MzU^^Ul)] (17) 

Using the MNS parametrisation in eq.(lO), the above three expressions (15) 
- (17) simplify to [16] 

= Y^^rCl2[c23Sl3Sl2^rl^31 " C23S 13 013^/^2^82 + t^rlt^T2Al] (18) 
^rC23Cl3[ci2C/rl^31 + 512^/^2^32] (19) 



dUe2 


1 


dt 


167r2 


dUeS 


1 


dt 


167r2 


dU^3 


1 


dt 


167r2 



^^13 ^ 1 , 2„ „2 

dt 167r2 

/irC23[-Sl2C^rl^l + Ci2C/r2^32] (20) 



^■^23 _ 1 ^2„2 

dt ~~ 167r2 
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2.3 Effect of scale-dependent vev 

In this section we modify the results of Ref.[16] by considering the running 
of the vev Vu{t) through the neutrino mass formula mi{t) — vl{t)Ki{t). This 
gives [15] 

d{lnmi) _ d{lnKi) ^ d{lnvu) 

dt ~ dt dt ^ ' 

where the RGE for in MSSM is given by 

^Vu = -^[—g'i + -gl-^h^]vu (22) 
dt 167r2^20^^ *^ " ^ ' 

The complete RGEs for neutrino mass eigenvalues are now given by [15] 

leading to the approximate solution, 

9 9 

miito) = mi{tR)exp{—Igi + -7g2)exp(-2C/^.j^) (24) 

Similarly, for Method A with the inclusion of scale-dependent vev, Rq in 
eq.(4) is now replaced by 

i?o = exp{^I,i + hg2) (25) 

Now the top-quark dependent term has been cancelled out and this will 
certainly affect the overall magnitude of the neutrino masses, but not the 
mixing angles. With this modification the magnitudes of neutrino masses 
tend to increase with the decrease of energy scale. The enhancement factor 
in the magnitude of neutrino masses through RGEs is calculated as 

/ = expi^h - - ^/,2) (26) 

which gives a positive numerical value greater than one even for large tan (5 
values (= 55). 
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3 Numerical analysis and results 



For a complete numerical analysis of the RGEs for both Methods A and B 
presented in the previous section, we follow here two steps (a) bottom-up 
running[7] and (b) top-down running[10]. In the first step (a), the running of 
the RGEs for the third family Yukawa couplings {hf , ht, hr) and three gauge 
couplings (91,92,93) in MSSM , is carried out from top-quark mass scale (to) 
at low energy end to high energy scale Mr where B — L symmetry breaks 
down[7,15]. In the present analysis we consider the high scale Mr — lO^^GeV 
and take the large tan/3 input value ( tan/3 = 55). For simphcity of the cal- 
culation, we assume here the supersymmctric breaking scale at the top-quark 
mass scale = lnmf[7]. We adopt the standard procedure to get the values 
of gauge couplings at top-quark mass scale from the experimental CERN- 
LEP measurements at Mz, using one-loop RGEs, assuming the existence of 
a one-hght Higgs doublet and five quark fiavours below rrit scale[15]. Simi- 
larly, the Yukawa couplings are also evaluated at top-quark mass scale using 
QCD-QED rescaling factors in the standard fashion [15]. Wc present here the 
values of the Yukawa couplings and gauge couplings at two scales to = 5.159 
and tu — 29.954 as follows: 



to 


= 5.159 


tR 


= 29.954 


tR 


= 29.954 


hf- 


= 1.000165276 


ht-- 


= 0.849373996 


It= 


:0. 109299302 


h 


= 0.866235097 




= 0.641425729 




=0.0759593919 




= 0.555977506 


hr- 


=0.679428339 


Ir = 


=0.0620016195 


91 


= 0.463751 


91= 


=0.626455069 


hi 


= 0.0544101298 


92 


= 0.6513289 


92 = 


=0.708234191 


h2- 


=0.0803824514 


93 


= 1.1891966 


93 = 


=0.784410238 


^3- 


= 0.141634345 



The values of the integrals are estimated between the two limits (to^i?)- 
In the second stage (b), the running of three neutrino masses (mi, 777.2,777.3) 
and mixing angles (si2, S13, S23) is carried out together with the running of 
the Yukawa and gauge cuoplings, from the high scale tu to low scale to- In 
this case wc use the values of Yukawa and gauge couplings evaluated earlier 
at the scale tu from the first stage runnung of RGEs. In principle one can 
evaluate neutrino masses and mixing angles at every point in the energy 
scale. 

We present the results of our numerical analysis in Tables 1-4. First we 
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check the stabihty of the neutrino mass models under radiative corrections. 
Table 1(a) and (b) give the values of neutrino masses and mixing angles at 
high and low scales for three neutrino mass models - hierarchical, inverted 
hierarchical and degenerate models (ses Appendix- A,B). For a check on nu- 
merical consistency, we evaluate the quantities for both Methods A and B 
outlined in the previous section. Both methods nearly give consistent results 
to about a factor of 13% discrepancy. Only the hierarchical model (Type 111) 
and the inverted hierarchical model (IIB) with opposite CP parity, are found 
to be stable under RG analysis in MSSM. In fact, both Amf2 and Am^g 
are slightly reduced, whereas the mixing angles are slightly increased with 
the decrease in energy scale. But there is no substantial change in both pa- 
rameters, which may spoil the good predictions already achieved at the high 
energy scale. In case of inverted hierarchical model with same CP parity ( 
IIA), the solar angle (^12) is not stable under radiative corrections. Similarly, 
in all the three degenerate models (I-A,B,C) both solar and atmospheric mass 
scales as well as solar and atmospheric angles, are not stable under radiative 
corrections. 

In Table 2 we analyse the MPR conjecture [12] which states that quark 
and lepton mixing angles are identical at high energy scale, and large solar 
and atmospheric neutrino mixing angles together with the small reactor an- 
gle, can be understood purely as a result of RG evolution provided the three 
neutrino masses are quasi-degenerate and have the same CP parity. In the 
present numerical analysis in Method B, both Yukawa and gauge couplings 
are running together with neutrino masses and mixings. It appears that in 
MPR paper [12] only neutrino masses and mixings are running. We present 
here three set of readings ( cases (i),(ii),(iii)) which give good radiative mag- 
nifications. However in case (iv) we use the same input values quoted in 
MPR paper but we observe less magnification contrary to MPR paper [12]. 
This discrepancy may be due to the running of gauge and Yukawa couplings 
along with neutrino masses and mixings in the form of coupled equations 
in the present analysis. In short, MPR conjecture[12] is well verified. The 
analysis presented here in Table 2 is the results of Method B, but we find 
that both methods (A and B) give consistent results. 

In Table 3 we supply some new results connected to two more conjectures 
on radiative corrections in MSSM. JM conjecture [11] specifies that radiative 
corrections can generate the neutrino mass-squared difference required for 
the large mixing angle (LMA) MSW solution to the solar neutrino problem 
if two of the three neutrino masses are assumed to be exactly degenerate 
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(m, —m, m') at high energy scale, and also if f/es at high scale is non-zero. 
We have shown in Table 3 that for a limited range of non-zero values of m', 
it is also possible to generate LMA MSW solution even if Ues is zero at high 
scale. If both m' and Ue3 are zeros at high scale, then it is not possible to 
generate LMA MSW solution. Case (ii) in Table 3 is interesting in the sense 
that it gives exception to JM conjecture by taking Ues = and m' — 0.07eV. 
At low energy we get Aml^ = 5.23 x lO^^eV^ and sin = 0.7042421 which 
is high but lies in the light side of the data (tan ^12 < 1) [21,22,23]. Fig.l 
shows the evolution of the three mixing angles, and the CHOOZ angle sin 6*13 
is generated through radiative corrections. Fig. 2 presents the evolution of 
the three neutrino masses, and the solar mass scale generated through ra- 
diative corrections at low scale, is then demonstrated in Fig. 3. It is still 
necessary to tone down the solar mixing angle through further fine tuning. 
As an example, using S12 = 0.7 at high scale, one can get S12 = 0.697 at lower 
scale. We also discuss MST conjecture [14] which states that starting from 
bimaximal mixings at high scale, radiative corrections due to the r- Yukawa 
coupling leads to solar angle towards the dark side at low scale, (tan ^12 > 1). 
The results presented in Table 3 also show that MST conjecture [14] is not 
always valid for the cases (i)-(iv) discussed here. 

Finally, as expected, the effect of running vev in neutrino mass formula 
leads to the increase of neutrino masses with the decrease of energy scale [15]. 
In Table 4 we calculate the neutrino masses and mixings with (I) or without 
(II) the effect of running the vev, at low energy scale. There is a factor of 
about 1.91 higher in case of low-energy neutrino mass values obtained with 
the running effect of vev. Figs. 4 and 5 show the evolution of neutrino masses 
for the above two cases (II) and (I) respectively. 
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Table 1(a): Running of neutrino masses and mixing angles from high 
scale Mr = lO^'^^GcV to top-quark mass scale rrit = 175GeV in MSSM for 
hierarchical (Type III) and inverted hierarchical (II) models {rriLL collected 
from Appendix B). Methods A and B are explained in the text. 



Type 


Item 




(A): fi = mt 


(B): fi = mt 


III 


mi 


0.00336eV 


0.0029358eV 


0.002502eV 




1712 


0.007357eV 


0.0063875eV 


0.00534eV 




ma 


0.057013eV 


0.0479200eV 


0.040322eV 




S23 


0.65630 


0.6632 


0.68441 




Sl3 


0.07358 


0.08141 


0.07911 




S12 


0.5838 


0.58807 


0.608521 




Aml^ 


4.28 X IQ-^eV"^ 


3.22 X IQ-^eV"^ 


2.23 X IQ-^eV^ 




Amis 


3.20 X IQ-^eV^ 


2.20 X IQ-^eV"^ 


1.60 X IQ-^eV'^ 


IIB 


mi 


-0.070445eV 


-0.0613552eV 


-0.05171696eV 




1712 


0.070977eV 


0.0618069eV 


0.0520974eV 






0.0005324eV 


0.0004482eV 


0.000378eV 




S23 


0.7071 


0.68486 


0.685182 




Sl3 


0.0 


0.00045 


0.000445 




Sl2 


0.7057745 


0.70581 


0.7057745 




Am?2 


7.52 X IQ-^eV^ 


5.56 X 


4.22 X 






1.96 X 10-^:V'- 


3.76 X 10-\:l'- 


2.67 X 10-'\:]"^ 


IIA 


mi 


0.0497257eV 


0.042245eV 


0.0356237eV 




m2 


0.0500693eV 


0.0446987eV 


0.037661eV 




m3 


0.000005eV 


0.000004eV 


0.000003eV 




S23 


0.707107 


0.68486 


0.68518 




Sl3 


0.0 


0.0 


0.0 




S12 


0.465 


0.99874 


0.998697 






3.43 X IQ-^eV^ 


21.13 X lO-'eV^ 


14.93 X IQ-^eV'^ 




Amis 


2.47 X 10-^6^2 


2.00 X lO-^eV^ 


1.42 X lO-^eV^ 
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Table 1(b): Running of neutrino masses and mixing angles from high 
scale Mr = 10^'^GeV to top-quark mass scale rrit = 175GeV in MSSM for 

degenerate models {ttill collected from Appendix B). Methods A and B are 
explained in the text. 

Type Item fj, — Mr [A): ji — rrit (B): /j, — rut 

~1A 0.396484eV 0.326001eV 0.275101eV 

m2 -0.396532eV -0.345072eV -0.290815eV 

ma 0.4eV 0.396941eV 0.3007399cV 

S23 0.707107 0.99885 0.998810 

si3 0.0 0.556 0.558 

si2 0.70931 0.853267 0.85202 

Amj^ 3.81 X IQ-^eV^ 12.80 x IQ-^eV^ 8.89 x IQ-'^eV^ 

Amja 2.76 x lO^^eV^ 3.85 x IQ-^eV^ 5.87 x IQ-^eV^ 

~1B 0.396841eV 0.315833eV 0.266586eV 

m2 0.396891eV 0.35615eV 0.3000754eV 

ms 0.4eV 0.357618eV 0.301312596eV 

S23 0.707107 0.99950 0.99949 Table 

si3 0.0 0.00854 0.011254 

si2 0.459701 0.99999 0.99999 

Am?2 3.97 X lO-^e^^ 27.09 x lO-^^e^^ 18.98 x lO-^^e^^ 

Amja 2.48 x lO'^eV^ 1.05 x IQ-^eV^ 0.74 x lO'^eV^ 

~lC rni 0.396841eV 0.334761eV 0.2823071cV 

1712 0.396891eV 0.35615eV 0.3000757eV 

ma -0.4eV -0.337397eV -0.284533eV 

S23 0.707107 0.7071 0.70722 

si3 0.0 0.0 0.0 

si2 0.4597 1 0.999999 

Am?2 3.97 X IQ-^eV^ 14.78 x IQ-^eV^ 10.35 x IQ-^eV^ 

AmL 2.48 X lO'^eV^ 13.01 x IQ-^eV^ 9.09 x lO^^eV^ 
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2: Analysis on MPR conjecture [12] related to radiative magnification on solar 
and atmospheric mixing angles at low scale. The parameters (m^ 2 3, 13 12) 
are defined at high scale Mr = lO^^GeV and others are defined at low scale 
rrit = 175GeV. Cases (i)-(iv) include difi^erent sets of arbitrary input param- 
eters. 



Case (i) Case (ii) Case (iii) Case (iv) 



m? = 0.5288446eV 


0.5288446eV 


0.5288446eV 


0.3682eV 


= 0.5309554eV 


0.5309554eV 


0.5309554eV 


0.370eV 


m° = 0.6eV 


0.6 eV 


0.6eV 


0.421eV 


s% = 0.0311 


0.03 


0.031 


0.038 


s% = 0.005 


0.0 


0.004 


0.0025 


s$2 = 0.22 


0.2 


0.215 


0.22 


mi = 0.399836034eV 


0.3997283eV 


0.39677897eV 


0.2783897eV 


m2 = 0.399989337eV 


0.4001846eV 


0.396988094eV 


0.2791149eV 


ma = 0.403035134eV 


0.40294793eV 


0.39996222eV 


0.28240252eV 


S23 = 0.709379 


0.71158975 


0.71025485 


0.467712 


si3 = 0.135225 


0.071834 


0.12237 


0.05389 


S12 = 0.656099 


0.6637173 


0.670588 


0.32458 


Am?2 = 12.26 X IQ-^eV^ 


36.50 X 10-^6^2 


16.60 X IQ-^eV^ 


40.43 X IQ-^eV^ 


Amig = 2.45 X 10~^eV^ 


2.22 X lO-W^ 


2.37 X 10-^6^2 


1.85 X IQ-^eV^ 
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Table 3: Analysis on JM conjecture [11] related to radiative magnification 
on solar mass scale at low scale. The parameters (m'j' 2 3, ^3 12) are defined 
at high scale M/j = lO^^GeV and others are defined at low scale rrit = 
175GeV. Cases (i)-(iv) include different sets of arbitrary input parameters. 



Case (i) Case (ii) Case (iii) Case (iv) 



ml = -0.095eV 


-0.095eV 


-0.08 eV 


0.08eV 


= 0.095eV 


0.095eV 


-0.08eV 


-0.08eV 


ml = O.OeV 


0.07eV 


0.04eV 


O.OeV 


= 0.707107 


0.707107 


0.707107 


0.707107 


s?3 = 0.1 


0.0 


0.1 


0.0 


s°2 = 0.707107 


0.707107 


0.707107 


0.707107 


mi = -0.0693038eV 


-0.06968284eV 


-0.0583183eV 


-0.058726eV 


m2 = 0.07013016eV 


0.070057414eV 


0.0591391eV 


0.058726eV 


ma = O.OeV 


0.049511038eV 


0.0283987eV 


O.OeV 


S23 = 0.6852 


0.638701 


0.6745 


0.6852 


si3 = 0.0969 


0.087023 


0.1329 


0.0 


S12 = 0.705008 


0.7042421 


0.70341 


0.707107 


Amj^ = 11.52 X IQ-^eV^ 


5.23 X IQ-^eV^ 


9.64 X IQ-^eV^ 


0.0e1/2 


Amia = 4.92 x lO'^eV^ 


2.46 X lO-^eV^ 


2.69 X lO-^ey^ 


3.45 X IQ-^eV^ 
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Table 4: Running of neutrino mass at low scale with (I) and without (II) 

the effect of scale- dependent vev w^. 



Item 


/i = MR 


/i = mt(II) 


jj, = mt{l) 


mi 


0.528844eV 


0.399836034eV 


0.762722015eV 


1712 


0.5309554eV 


0.399989337eV 


0.763014257eV 




0.6eV 


0.403035 134eV 


0.768824458eV 




0.0311 


0.709379 


0.709379 


Sl3 


0.005 


0.135225 


0.135225 




0.22 


0.656099 


0.656099 


Am?2 


2.24 X IQ-^eV^ 


12.26 X 10-^6^2 


44.588 X 10-^6^2 


Amig 


7.81 X 10-2eT/2 


2.45 X lO-^eV^ 


8.90 X lO-^eV^ 



4 Summary and Discussion 

We summarise the main points in this work. First we briefly review the main 
points of the formalism based on two approaches on the evolution of RGEs 
of neutrino masses and mixings. The first one (A) deals with the running of 
the whole neutrino mass matrix from which one can extract mass eigenvalues 
and mixings at any particular energy scale, whereas in the second approach 
(B) the three neutrino mass eigenvalues and the three mixing angles are run- 
ning directly. Detailed numerical analysis shows that both approaches agree 
up to a discrepancy of a factor of 13% in mass eigenvalues. The predictions 
on mixing angles are almost consistent in these two approaches. Using both 
approaches we show that hierarchical model (III) and inverted hierarchical 
model with opposite CP parity (IIB) are stable under radiative corrections. 
The evolution of sin 612 is very fast in the inverted hierarchical model with 
same CP parity (HA), and hence the model is not stable. We also verify 
the MPR conjecture [12] in which radiative magnification of solar and atmo- 
spheric mixings are possible in case of nearly degenerate model with same 
CP parity. We find that runnings of masses and mixings with and without 
the running of gauge and Yukawa couplings, give reasonably different magni- 
fications. However such radiative amplification generally involves a delicate 
fine-tuning of the initial conditions which are to some extent unnatural [17]. 
However such problems are not there in the theory of neutrino masses de- 
rived from Kahler potential in supersymmetric model, and neutrino mixings 
angles can easily be driven to large values at low energy as they approach 
infrared pseudo-fixed points at large mixing [17]. 



14 



We also study JM conjecture [11] which specifies the radiative generation 
of solar scale in exactly two-fold degenerate model having opposite CP parity 
(m, — m, ms) and non-zero values of f/es- Our numerical analysis shows that 
the same radiative generation of solar mass scale is also possible with the 
conditions C/es = 0.0 and non-zero value of 1113. We also discuss the MST 
conjecture[14] which states that starting from bimaximal mixings at high 
scale, radiative corrections lead to the solar angle towards the dark side of 
the data at low energy scale. We show that MST conjecture is not always 
valid. We suggest further generalisation of JM conjecture and this will be 
reported in subsequent communication [25]. Finally, the effect of running the 
vev in neutrino mass renormalisation is discussed and it is observed that 
neutrino mass increases with the decrease of energy scale when we include 
the running of vev. This gives magnification of a factor of 1.91 in neutrino 
masses at low energy scale, compared to the values calculated without the 
running vev effect. Numerical analysis in three generations with arbitrary 
CP violating phases is of great interest [18] and there has been a possibihty 
that within a restricted range of the physical parameters including phases, 
the degenerate models are found to be stable under radiative corrections in 
MSSM [26]. As emphasised before, oversimplifications of analytic expressions 
and a departure from the simultaneouse running of Yukawa and gauge cou- 
plings along with neutrino masses and mixings, may have a danger of getting 
misleading conclusions. 



Acknowledgements 

N.N.S. thanks the High Energy Physics Group, International Centre for The- 
oretical Physics, Trieste, Italy, for kind hospitality during the course of the 
work. 



15 



Appendix A 



The zeroth-order left-handed Majorana neutrino mass matrices with tex- 
ture zeros, niLL, corresponding to three models of neutrinos, viz., degenerate 
(Type [I]), inverted hierarchical (Type [II]) and normal hierarchical (Type 
[III]). These mass matrices are competible with the LMA MSW solution as 
well as maximal atmospheric mixings [19]. 

ri-i diaa 

Type rriLL 



[lA] 



/ 



1 1 _ 1 

V2 2 2 

J_ _ 1 1 

\ V2 2 2 / 



mo Diag{l, —1, l)mo 



[IB] 



1 
1 I mo 
1 



Diag{l, 1, l)mo 



[IC] 



1 
1 
1 



mo 



Diag{l, 1, — l)mo 



[IIA] 



/ 1 
I I I mo 



Diag{l, 1, 0)mo 



[IIB] 



/ 1 1 \ 

1 

VI oj 



mo Diag{l, —l,0)mQ 



[III] 



/ 




\ 

1 _ 1 

V -i - 

\ " 2 2 



mo Diag{0, 0, l)mo 
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Appendix B 

The left-handed Majorana mass matrix ttill for three different models of 
neutrinos presented in Appendix A. These results are collected from Ref.[19]. 

Degenerate model( Type [lA]): 



rriLL = 



/ (-25i + 282) 



0.4 



Degenerate model (Type [IB]): 



(1 - 25i - 282) -Si 

rriLL = I -5i (1 - -S2 I 0.4 

-5i -52 (1 

Degenerate model (Type [IC]): 




rriLL 



/ (1 - 2Si - 2S2) -Si 

-Si 

V -^1 



Invereted hierarchical model(Type [HA]): 



-Si 



-S2 (1 - ^2) I 0.4 
1 - ^2) -S2 



[1 - 2e) -e -e 



rriLL 



— e 



\ 0.05 



Inverted hierarchical model(Type [HB]): 





( ° 


1 


1 


■niLL = 


1 









I 1 





A3 



Hierarchical model (Type [III]: 
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-A^ A A3 \ 
A 1-A -1 0.03 
A^ -1 1-AV 



The values of the parameters used are: Type lA: 6i = 0.0061875, S2 — 
0.0030625, Type [IB] and [IC]: di = 3.6 x 10-^ 62 = 3.9 x 10-^ Type 
[IIA]: 77 = 0.0001, e = 0.002, and Type [JIB] and [III]: A = 0.22. All neutrino 
masses are in eV. 
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Fig.l Evolution of the three mixing angles with energy scale in JM conjec- 
ture [11]. sin 023, sin6'i3 and sin 12 are represented by solid line, dashed- 
line and dotted-line respectively. 
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Fig. 2 Evolution of the three neutrino mass eigenvahies in JM conjec- 
ture [11]. mi, 1712 and ma are represented by dashed-hne, solid- line and 
dotted-line respectively. 
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5 10 15 20 25 30 

t = In/t 

Fig. 4 Evolution of the three nearly degenerate neutrino masses in MPR 
conjecture[12] (without the effect of scale-dependent vev in II). ma, m2 and 
mi are represented by solid line, dashed-line and dotted-line respectively. 
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5 10 15 20 25 30 



Fig. 5 Evolution of the three nearly degenerate neutrino masses in MPR 
conjecture[12] (with the effect of scale-dependent vev in I), m^, m2 and 
mi are represented by solid line, dashed-line and dotted-line respectively. 
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